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ABSTRACT 

NGC 4261 (3C 270) is a low-luminosity radio galaxy with two symmetric 
kiloparsec-scale jets. Earlier Hubble Space Telescope observations indicated 
the presence of a hundred-parsec scale disk of cool dust and gas surrounding 
a central, supermassive (~ 4.9 x 10 8 M & ) black hole. The recent detection 
of free-free radio absorption by a small, geometrically-thin disk, combined 
with earlier studies of the disk's large scale properties, provide the strictest 
constraints to date on the nature of the accretion process in this system. We 
show here that a supersonic disk, illuminated by the active galactic nucleus 
(AGN), can not only account for the observed radio shadowing, but can also 
produce the optical broad lines emitted from this region. At large radii, the gas 
is optically-thin because the ram pressure due to turbulence is much larger than 
the thermal pressure of the gas. At smaller radii, but beyond a critical radius 
r c , line cooling dominates over gravitational dissipation and the gas is effectively 
cooled down to temperatures below 10 4 K. Within r c , however, heating due to 
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the release of gravitational energy overwhelms line cooling and the gas falls onto 
the unstable portion of the cooling curve. Because cooling is quite inefficient 
under these conditions, the plasma is heated very quickly to a temperature close 
to its virial value as it falls toward the central engine. Thermal pressure of the 
gas dominates the turbulent ram pressure at a radius ~ (2/3) r c , below which 
the flow probably becomes advection dominated. The disk is optically-thin to 
UV and X-ray radiation within r c , so the ionizing radiation from the AGN is 
preferentially absorbed near r c , affecting the disk structure significantly. To 
include the ensuing photoionization effect, we have used the algorithm Cloudy 
with additional heating introduced by gravitational dissipation to calculate the 
temperature profile and line emission from the disk in a self-consistent manner. 
The results of our model calculation are consistent with current multiwavelength 
observations of the disk in this source. 

Subject headings: accretion — black hole physics — Galaxy: NGC 4261— 
hydrodynamics — turbulence — line: formation 

1. Introduction 

The development of supermassive black hole accretion theories has always been limited 
by observations (Pringle 1981). This has been changing quickly during the past few years, 
with increasingly higher resolution observations being made over a wide frequency range. 
For example, the Chandra X-ray Observatory {Chandra) and the Hubble Space Telescope 
(HST) both have sub-arcsec scale resolutions, while the Very Large Baseline Interferometer 
( VLB!) can resolve structures on sub-milliarcsec scale. By combining multiwavelength 
observations from these different instruments, we can obtain very strict constraints on the 
nature of accretion around supermassive black holes. 

While the innermost structures of AGNs remain unresolved (however, cf. Melia & 
Falcke 2001 regarding the Galactic center), the current data are good enough to constrain 
the properties of the large-scale accretion disks that are presumably feeding the central 
engines. A better understanding of these large-scale disks is not only tenable, but will also 
provide significant implications on the nature of the AGNs. Most of the relevant theoretical 
work to date is based on the standard a-disk prescription (Shakura & Sunyaev 1973), in 
which the accretion disk at large radii is supported by the thermal pressure of the gas (see, 
e.g., Hubeny 1990; Hure et al. 1994; Collin & Hure 1999). This may not be the entire story, 
however, as some observations suggest that turbulence is prominent in the disks of many 
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radio galaxies (Ferrarese, Ford & Jaffe 1996; Ferrarese & Ford 1999; van Langevelde et 
al. 2000). In addition, the inner disk is irradiated (Maloney, Hollenbach & Tielens 1996), 
which should affect, e.g., its physical properties and line emission characteristics. 

At a distance of only 30 Mpc, the FR I radio galaxy NGC 4261 provides us an 
exceptional opportunity to investigate these new effects. Harboring a 4.9 x 10 s M Q black 
hole at its center (Ferrarese et al. 1996), NGC 4261 is well known for its kiloparsec-scale 
jets (Birkinshaw & Davies 1985) and its hundred parsec-scale dusty disk imaged by HST 
(Jaffe et al. 1993). This disk possesses an atomic surface density of ~ 5 x 10 20 cm~ 2 , as 
inferred from its optical depth assuming a Galactic dust/gas ratio. On a smaller scale, 
a thin circumnuclear disk of radius ~ 6 pc was revealed by radio line absorption studies 
(Jaffe & McNamara 1994; van Langevelde et al. 2000). This portion of the disk has an H I 
column density of ~ 2.5 x 10 19 T sp cm" 2 , where T sp ~ 100 K is the spin temperature of 
the gas. At an even smaller scale, free-free absorption of radio emission from the core by 
a geometrically-thin disk of radius ~ 0.3 pc was detected (Jones et al. 1997, 2000, 2001; 
Piner, Jones & Wehrle 2001). This recent finding adds significant content to the accretion 
picture because this radius corresponds to only 6400 r^, where the Schwarzschild radius 
r s = 2GM/c 2 ~ 1.45 x 10 14 cm for the central black hole in NGC 4261. 

Although these observations of the accretion disk extend to sub-parsec scales, a 
fully self-consistent accretion model does not yet exist, particularly for the small free- free 
absorbing disk (Sudou et al. 2000). Jones et al. (2000) did have success in attributing 
the observed free-free absorption to a standard a-disk. Unfortunately, to match the 
observations, the model requires a dissipative stress that is much greater than that provided 
by subsonic turbulence. A viable alternative for (or addition to) this stress must be found 
in order for the model to be fully self-consistent. 

We suggest here that supersonic turbulence is the responsible agent. Observations of 
NGC 4261 provide direct evidence that turbulence is prominent in the disk from 100 pc to 
about 6 pc, with the inner boundary set by observational limitations. At 100 pc, the width 
of narrow lines from the disk suggest a turbulent velocity of 70 km/s (Jaffe et al. 1996), 
while at 6 pc, the turbulent velocity estimated from the H I absorption profile is about 
130 km/s (Jaffe & McNamara 1994; van Langevelde et al. 2000). Strong turbulence in 
the accretion disk is understandable given that the disk is embedded within the core of 
the galaxy, where it is subject to disruption by strong stellar winds. The dissipation of 
turbulence is regulated by the cascade in energy from large scales to small scales, so the 
total energy associated with the turbulent motion decays roughly as a power-law in time 
(Lazarian 2003), even in the case of supersonic turbulence (Smith, Mac Low & Zuev 2000a). 
This slow turbulence decay suggests that the supersonic structure might persist even to 
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smaller radii. 

In such a supersonic disk, the accreting gas is optically-thin and gravitational energy- 
dissipation will overwhelm line cooling at a critical radius r c , within which the gas is 
effectively heated to a temperature close to its virial value (Pringle 1981). Thermal pressure 
of the gas will then dominate the turbulent pressure at r ~ (2/3)r c , so the supersonic 
flow appears to terminate at this radius and the gas may make a transition into an 
advection-dominated accretion flow thereafter (Narayan & Yi 1994). On larger scales, the 
supersonic structure can be inferred directly from observations. It is interesting to note 
that with a stress parameter a ~ 1, the value of r c is roughly 0.3 pc, which corresponds 
to the inferred size of the free-free absorbing region. Thus the observed shadowing of the 
nascent radio jet may be produced by the ionized gas near r c . 

We further note that newly released Chandra data also hint that the nucleus of 
NGC 4261 may be weakly-active (Peter, Zezas & Fabbiano 2001), consistent with earlier 
ASCA and R OSA T observations (Matsumoto et al. 2001; Worrall & Birkinshaw 1994). The 
corresponding ionization parameter at 0.3 pc (assuming that the radiation is isotropic) is 
S 5, where S = F ion / n to tkTc, F ion is the ionizing flux between 1 and 10 3 Ryd, and 
n to t is the total number density including all species. Note that S roughly corresponds to 
the ionizing-photon-to-gas pressure ratio. A value of S ^ 1 indicates that the gas may be 
effectively photoionized. This photoionization will significantly affect the profiles of the 
lines emitted by the supersonic disk. High-resolution HST imaging of the central ionized 
gas has detected broad emission lines from a region smaller than 7 pc in radius (Ferrarese 
et al. 1996). In this paper, we show that a supersonic disk, irradiated by the AGN, can 
not only account for the observed radio absorption feature, but also reproduce these broad 
emission lines. We also show that its large-scale extension is consistent with the radio line 
absorption observations. 

The physical conditions present in the broad line regions (BLRs) of this class of radio 
galaxy are likely to be very different than those present in the more-luminous Seyfert 
galaxies and quasars. In both scenarios, the line emission is produced by gas having 
temperature on the order of a few times 10 4 K, with line broadening caused by bulk motions 
of the emitting plasma. However, while the BLR in Seyferts and quasars is primarily 
photoionized, our model predicts that heating due to gravitational dissipation plays a 
significant role in the low-luminosity radio galaxies. 

Our modeling also indicates a very different BLR geometry for these galaxies, with 
the line emission coming directly from the turbulent accretion disk. While accretion 
disks may play an important role in some BLR models of Seyferts and quasars, the line 
emission is not believed to originate in the disk itself, instead being associated with, e.g., 
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magnet ohydro dynamic flux tubes (Emmering, Blandford, and Shlosman 1992) or radiation- 
and gas-pressure driven winds rising from the surface of the disk (Murray et al. 1995). The 
accretion disk does not even play a direct role in many successful BLR models for these 
objects, e.g., the bloated star (Alexander & Netzer 1994) or accretion- shock cooling models 
(Fromerth & Melia 2001). 

It is also noted that the density of the BLR gas in these galaxies is much lower than in 
traditional AGNs. Because broad forbidden emission lines are absent in Seyfert and quasar 
spectra, we know that the density of their line-emitting gas is n ^ 10 7 cm™ 3 . However, the 
spectra of some low-luminosity radio galaxies contain strong broad forbidden emission lines 
(e.g., [N II], [S II]), indicating gas densities n ^ 10 5 cm™ 3 (see, e.g., Ferrarese & Ford 1999; 
Barth et al. 1999). The lower density BLR gas may be directly tied to their lower mass 
accretion rates and resulting luminosities. 

We first discuss the general properties of a supersonic disk in § 2. Its consistency with 
the multi-wavelength observations of NGC 4261 is discussed in § 3. Detailed modeling of 
the line emission from this supersonic disk, irradiated by its AGN, is given in § 4. In § 5, 
we discuss the implication of this work on the nature of AGNs in radio galaxies. 



Because it is thought that the large-scale dusty disk feeds the inner accretion of the 
AGN (Jaffe et al. 1996), the detection of a prominent gap in the radio core of NGC 4261 
prompted Jones et al. (2000, 2001) to introduce a free-free absorbing disk to account for the 
observed features. To produce the radio absorption, the inner disk must have an emission 
measure of nil ~ 3 x 10 8 T^ 2 pc cm™ 6 , where n e is the electron density, / is the path 
length through the absorbing gas, and T 4 is the electron temperature in units of 10 4 K. If 
dominated by thermal gas pressure, however, the optically thin a-disk proposed by these 
authors cannot account for this radio shadowing in a self-consistent manner. In particular, 
we show below that the implied stress parameter a is too large for the flow to be subsonic. 

For a gas-pressure-dominated accretion disk, angular momentum conservation requires 
that v ~ — (2/3) f r r at large radii, where v = aH(r) c s (r) is the kinetic viscosity, a is the 
stress parameter (assumed to be radius- independent and less than one), v r is the radial 
velocity of the flow, H(r) is the scale height of the disk, and c s (r) is the local sound speed. 
From the expressions given in Jones et al. (2000), one can obtain: 



2. A Supersonic Accretion Disk Model 
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where v k = (GM/r) 1 / 2 is the Keplerian velocity, M 8 is the black hole mass in units of 
1O 8 M , M_ 3 = M/(10 -3 Mq/jt) gives the accretion rate, and r 18 ~ 1 is the distance of the 
radio absorbing gas from the black hole in units of 10 18 cm. In fact, it can be shown that 
the disk cannot be supported by thermal pressure of the gas alone since the corresponding 
accretion rate is: 



M ~ 7.4 x 10 19 a 



2.0 x 10 5 cm- 3 / V10 4 K 



/ T \ 3 / 2 

/ J o \ 



X (4.9 x 10 8 M ) (l0 18 °cm) gS ' ' (2) 

where quantities with subscript o denote the corresponding values inferred from the radio 
absorption features, and we have approximated the line-of-sight path length through the 
disk with 2#/sin(26°) (Jones et al. 2001), where H = (R g Tr 3 / fiGM) 1 ^ 2 , R g is the gas 
constant and /i = 0.5 is the molecular weight. Such an accretion rate is so small that the 
implied emission efficiency for converting accreted rest-mass energy into X-ray emission 
from the nucleus of NGC 4261 is about 3.0 /a (Worrall & Birkinshaw 1994; Matsumoto 
et al. 2001), which is unacceptable for any a < 1 in a thermal pressure supported disk. 
Sudou et al. (2000) further discussed an optically-thick standard disk and an optically-thin 
disk cooled by free-free emission. Neither of these alternatives can account for the observed 
shadowing. 

Meanwhile, several observations indicate that turbulence dominates the gas motion 
at large scales (Jaffe et al. 1996; van Langevelde et al. 2000), where supersonic turbulence 
may be produced by stellar winds in the core of the galaxy (Jaffe et al. 1996). As the 
gas falls in toward the black hole, however, the strength of the residual turbulence in the 
disk is determined by balancing shock dissipation with the turbulence enhancement due 
to gravitational energy dissipation. In a fully developed turbulent flow, energy is mostly 
associated with large-scale gas motion. The decay rate of the total turbulent energy is then 
determined by the rate at which energy cascades from large scales to small scales. The 
turbulent energy density per unit mass, e, is therefore given by 

i = e/r, (3) 

where the dot denotes a derivative with respective to time and r = Lje 1 ! 2 gives the 
characteristic time scale for transporting energy from large to small scales. The characteristic 
length scale, L, of the system depends on the turbulence generation mechanism. Assuming 
L is independent of time, we have e oc t~ 2 . In their numerical simulations of the decay of 
supersonic turbulence, Smith et al. (2000a) found that the total energy of the turbulence 
decays roughly with a power law in time, with an index —1.5. 
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We point out that as the large-scale supersonic gas flows inward, the dissipative stress 
is provided primarily by its turbulent motion. So the dissipated gravitational energy will 
go first into producing turbulence. This turbulence generation mechanism can, in principle, 
keep the disk supersonic even on smaller scales (Smith, Mac Low & Heitsch 2000b). A 
complete model addressing the evolution of this supersonic accretion flow, however, would 
require a detailed numerical simulation beyond the scope of this paper. 

In the following discussion, we will assume that the energy associated with the 
turbulent motion is in sub-equipartition with, and proportional to, the gravitational energy 
of the gas. We therefore have 

v 2 t = fv 2 k , (4) 

where / is a "turbulence parameter" , which can be fixed by the observed disk properties at 
large radii, and v t is the turbulent velocity in the flow. The supersonic disk model (in which 
the stress is dominated by turbulence) then gives: 

v = a H v t , (5) 

H = fV*r, (6) 
v 

v r = —1.5 — 

r 

= -1.5 afv k , (7) 
M = 6vr a f /2 r 2 v k nm p . (8) 

In Eq. (8), n(r) is the baryon number density averaged in the vertical direction of the 
accretion disk and m v is the proton mass. Because the H I absorption feature indicates that 
Vt — 0.21 Vk at r ~ 6 pc from the nucleus, we obtain / ~ 0.045. The corresponding disk 
opening angle is 2 tan -1 / 1//2 ~ 24°. Then from Equations (6) and (8), we have 

H = 0.21 r , (9) 

/r \ 3 / 2 

n = n ^j . (10) 

It is interesting to note that the viscous time scale in such a disk (t V i S — r/v r by definition) 
is longer than the corresponding dynamical time scale (td = 2-wr /vk) by a factor of 2.4/a. 
Adopting the power-law decay rate from the numerical simulations (Smith et al. 2000a), 
one can show that the ratio of the turbulent energy to the dissipated gravitational energy 
is about 0.27a 3 / 2 . This is consistent with / ps 0.045 for a less than one. 

One should emphasize that only a small fraction (~ 4.5%) of the dissipated gravitational 
energy is left to drive the turbulent motion. Most of the dissipated gravitational energy 
is thermalized via shocks and effectively radiated away. Since we are mostly interested in 
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regions where the gas is partially ionized, this cooling is dominated by line emission. We will 
see that there is a critical radius r c in the disk, where the gravitational energy dissipation 
rate equals the line cooling rate. The gravitational energy dissipation rate is given by 



r = 



3 GMM 
8tt Hr 3 



From Equation (8), we have 
M = 
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Using this equation to eliminate M in Equation (11), we get 
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The line cooling rate of an optically-thin gas is given by 



A = n 2 A 



9.0 x 10~ 14 A_ 
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3 x 10 4 cm" 3 
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10 18 cm 



10 18 cm 



ergs cm 3 s 1 , 



where A is the cooling function (Krolik 1999) and A_22 = A/(10 22 ergs cm 3 s 1 ). 
In terms of the critical radius r c , we have 
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(12) 



(13) 



(14) 



(15) 



To account for the radio shadowing, the electron number density should be ~ 3 x 10 4 cm -3 
at ~ 0.3 pc, where we have assumed an inclination angle % = 64° for the disk (Jones et 
al. 2001). Thus, for a stress parameter a ~ 1, Equation (15) shows that the critical point is 
located within the radio shadowing region. It is natural to expect that the observed radio 
absorption should be produced by the fully ionized gas near r c , and we will see below that 
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ionizing photons from the central engine play a crucial role in determining the ionization 
state of the gas at this radius. 

HST observations indicate that the gas closest to the nucleus may be ionized by 
radiation from the central engine (Ferrarese et al. 1996). Further evidence for this can 
be shown by estimating the ionization parameter at r c . Recent high-resolution Chandra 
observations have uncovered a high-energy source at the nucleus of NGC 4261 with a hard 
band X-ray luminosity of L ps 7 x 10 40 ergs/s (Peter et al. 2001), confirming earlier ASCA 
results (Matsumoto et al. 2001). Meanwhile, the soft X-ray luminosity is ~ 10 41 ergs/s 
(Worrall & Birkinshaw 1994). So, assuming that the radiation is isotropic and that the 
gas temperature is ~ 10 4 K, we can estimate the value of the ionization parameter to be 
S ~ 5 for the radio- absorbing gas at ~ 10 18 cm (Krolik 1999). The radio-absorbing gas is 
therefore effectively photoionized. 

Although the radiation pressure dominates the thermal pressure of the gas in the radio 
absorbing region, the turbulent pressure there is 0.5 n m p v^ ~ 7.4 x 10~ 5 ergs cm -3 , 
which is much larger than the radiation pressure. So the disk structure, including the scale 
height and number density, is still determined by the turbulent motion. However, the gas 
temperature profile depends on the interaction between the disk and the radiation field. In 
the following sections, we will incorporate this photoionization effect using the algorithm 
Cloudy (Ferland 1996; Cloudy 96-/33) and show that the observed broad-line spectrum may 
also be fitted with this disk model in a self-consistent fashion. 



3. Line Emission from the Irradiated Supersonic Disk 

To model line emission from the irradiated disk, we first need to know the intensity 
and spectral energy distribution (SED) of the ionizing radiation from the central source. 
These are highly uncertain due to a combination of factors. First, there are uncertainties 
associated with intrinsic absorption and galactic contamination which limit the precision 
to which we can quantify the spectrum in the optical and X-ray bands. Second, the 
spatial distribution of photons is unknown — it is unclear whether the central engine 
emits isotropically; the disk at r c may "see" a different spectrum from what we observe. 
Third, and most importantly, we do not have observations (or even good estimates) of the 
continuum at UV energies. Generally, the shape and intensity of the UV continuum most 
strongly determine the properties of a photoionized gas. 

The correlations between radio and optical spectra have suggested a non-thermal 
nuclear source for FR I radio galaxies (Chiaberge, Capetti & Celotti 1999). Further 
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evidence for a power-law nuclear source in NGC 4261 comes from ROSAT and ASCA 
observations in the X-ray band (Worrall & Birkinshaw 1994; Matsumoto et al. 2001), which 
are confirmed by recent high spatial resolution Chandra observations (Rinn & Sambruna 
2001; Peter et al. 2001). In the following discussion, we will adopt a power-law spectrum 
in the UV-range to bridge the gap between the optical and X-ray data. The continuum is 
assumed to be emitted isotropically from a central, point-like source, and the intensity of 
the ionizing luminosity remains a free parameter in our model to reflect our uncertainties. 

We are mostly interested in the region at ~ 0.3 pc. However, to determine the 
spectrum of the incident radiation, we also need to know the gas distribution in the inner 
region since radiation from the core may be reprocessed by intermediate gases. Fortunately, 
according to the supersonic disk model, the disk will be heated up very quickly within r c , 
which makes the gas there optically-thin to UV and X-ray emission. One can show this by 
solving the energy conservation equation in the transition region between r c and the radius 
r t , where the thermal pressure of the gas equals the turbulent pressure and the supersonic 
accretion flow effectively terminates. 

For a supersonic accretion disk, we have the energy conservation equation: 



d_ 

dr 



1.5(1- f)Mv 2 k - 



5Mk B T 



-Air r H A . 
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The first term on the left-hand side corresponds to the gravitational energy dissipation and 
turbulent advection, while the second term is associated with the thermal energy of the gas. 
The right-hand side gives the cooling. The quantity A/n 2 = A is about 10~ 22 ergs cm 3 s _1 
near the peak (at ~ 10 4 K) of the cooling function (Krolik 1999). Assuming that A does 
not depend on radius, we have the following temperature profile in the transition zone: 
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In the second expression above, we have kept only the term proportional to [(r c — r)/r] 2 
and have neglected the high order terms in (r c — r)/r. Setting 2 /c B T(r t ) = / m v v\, we 
obtain the termination radius r t of the supersonic flow: 



(17) 
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1 + a/2 + (a + a 2 /4)V2 ' 



(18) 
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where 

a = oJfhn * 015 (oL) ■ (19) 

So r t ~ (2/3)r c , and T(r t ) ~ 3 x 10 7 K. Then the error introduced by neglecting the high 
order terms in Equation (17) is ~ 30%. In general, the termination radius r t is determined 
by solving the equation 

h, fro = 1 + f'\ (20) 



\rj 3(l-/)r, 

So / must be smaller than 3/8 to make the thermal pressure of the gas dominate at 
small radii. A larger value of / will cause the turbulent advection to cancel most of the 
gravitational energy dissipation and keep the disk supersonic. In addition, a more accurate 
expression for A in the transition zone will make the gas heated even faster because it is 
on the unstable portion of the cooling curve (Liu, Fromerth & Melia 2002). Within r t , the 
accretion pattern may transfer into an advection dominated flow (Narayan & Yi 1994). 

Within r c , the disk is optically-thin to UV and X-ray radiation, and because of the 
sharp increase in the gas temperature in this region, free-free absorption is also negligible. 
Due to the effects of photoionization by the AGN, however, the gas will be fully ionized well 
beyond r c . To incorporate this effect self-consistently, we will need to use the algorithm 
Cloudy with additional heating introduced by gravitational dissipation to calculate the 
temperature and the ionization state of the gas everywhere in the disk. 

Before calculating the line emission from the supersonic disk numerically, it may be 
illustrative to estimate the physical conditions of the broad-line region. Observations 
with HST have identified an ionized gas concentration in a spatially-resolved region 
with a FWHM of 17 pc. The broad emission lines are associated with the central 
O'.'l ~ 14 pc aperture position (Farraese et al. 1996). For the Ha line, the observed flux is 
3.11 x 10~ 15 ergs cm~ 2 s' 1 with a FWHM of 2500 km/s. From the electron number density 
inferred above and the Ha volume emissivity of ~ 2 x 10~ 25 ergs cm 3 s _1 (Osterbrock 1974), 
we can estimate the size of the emission region. In our model, most of the lines should be 
produced beyond r . The radius r of this region should therefore satisfy the condition 

4.8 x 10 38 ( -L.) V2 ( ^ ) 2 ( —^—) 3 In (r/r ) = 3.3 x 10 38 , (21) 

V0.045/ V3xl0 4 cm- 3 / V10 18 cm/ K / ' v ' 

implying that r m 0.6 pc. It is interesting to note that the Keplerian velocity at such a 
radius is 1900 km/s, which is consistent with the FWHM of the Ha line if one takes into 
account the turbulent motion in the flow and the inclination of the disk. Because the 
contributions to Ha from the various rings at different radii are comparable to each other 
between ~ 0.3 pc and ~ 0.6 pc, the double peaks that are usually produced by Keplerian 
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motion are here filled in, consistent with the observed spectrum. So a supersonic accretion 
disk with the properties described above appears to account for the broad emission lines 
rather naturally. 

4. Modeling the Broad Emission-Line Spectrum 

Several observational constraints already limit the parameter space available for fitting 
the broad emission lines. Images from HST show that the large-scale dusty disk has an 
inclination angle of 64° (Jaffe et al. 1993). On the other hand, the narrow- line kinematic 
data from the central region suggests an inclination angle of 69° for the ten-parsec scale 
disk. Both of these values are consistent with the viewed angle of the radio jet, which is 
~ 63° (Piner et al. 2001). We therefore adopt a value of i = 64° for the inclination angle of 
the supersonic disk. 

As stated in Section 2, the radio-line absorption observations fix the turbulent 
parameter at / = 0.045. To account for the radio shadowing near 0.3 pc, the electron 
number density should be around 30, 000 cm -3 . Because we expect hydrogen to be 
fully ionized in this region, we will choose n Q = 30, 000 cm~ 3 as the fiducial value for 
the scale-height-averaged number density of the disk at a radius r Q = 10 18 cm from the 
supermassive black hole. The remaining free parameters in the model are the ionizing 
luminosity L ion , and the stress parameter a, which determines the accretion rate and critical 
radius via Equations (12) and (15). 

4.1. Methodology 

We use Cloudy to calculate the line emissivity as a function of radius for r > r c , 
including the additional heating (Eq. 13) and varying the density with radius (Eq. 10) 
as predicted by our model. The line profiles are then generated by integrating over the 
appropriate range in radius. We accomplish this task by partitioning the disk into a large 
set of concentric tori, modeling the line emission from each torus as a set of N discreet 
clouds, with the understanding that we approach the continuous limit with very large N. 
The emission line flux from the population of N clouds in torus i as a function of wavelength 
is: 

N 

FfaX) oc 2nr i H{r i )Ar i ]T efo, Vj ) Afo, Vj , A) , (22) 
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where the factor 2n r\ H{r-i) Arj is the volume of torus i, e(r J -,v : ,) is the line emissivity 
including relativistic effects, and A(i\,-, Vj, A) is a convolution term relating a cloud's location 
and velocity to the observed line wavelength. The composite line flux from the entire disk 
is then obtained by summing the contributions from each torus. 

Taking the line-of-sight velocity to be i>i os , the effect of Doppler boosting is given by 
(Corbin 1997): 

f(r ' v) = eo(r) Xtt^jt) • (23 > 

where eo( r ) is the rest-frame emissivity (generated by Cloudy) and 



vL = t^ + c^l-^-l) (24) 

is the line-of-sight velocity including general-relativistic corrections. 

The observed wavelength A of line emission from a given cloud is related to rest-frame 
line wavelength Aq by 



A = A j 1 - I^f V2 (A(L _) (1+ :). (25) 



1 - Vlos/C 

where the correction terms on the right are the gravitational, Doppler, and cosmological 
redshifts, respectively. This relation, appearing as a delta- function A(i\,-,Vj,A) in Eq. (22), 
is used to map each cloud's emission into predetermined wavelength bins. 

The axis of the disk is taken to lie along the z-axis. Cloud positions and velocities 
are then determined via Monte Carlo sampling. Each cloud is given a random azimuthal 
position < < 2ir in the disk, and the velocity of each cloud is assigned according to 

Vj = v k + v t , (26) 

where v k = (GM/r) 1 / 2 is the Keplerian velocity, is the azimuthal direction vector (note 
that this depends on the value of 0), and \ t = v7 v k in a randomly-oriented direction. The 
line-of-sight velocity is then equal to v\ os = Vj • ni os , where ni os is the direction vector along 
the line of sight. 



4.2. Results of the Line Modeling 



We have modeled the continuum-subtracted Ha A6563 + [N II] AA6548, 6584 and 
[S II] AA6717, 6731 broad emission-line spectrum of Ferrarese et al. (1996). The parameters 
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L ion and a were varied to determine the best \ 2 fit to their nuclear pixel data, corresponding 
to a FWHM spatial resolution of O'.'l = 14 pc. Figure 1 shows the resulting best fit 
(xl = 11; v = 250 degrees of freedom), having parameters L ion = 4 x 10 40 erg/s and 
a = 0.4. The model nicely reproduces the broad wings and steep core of the Ha + 
[N II] composite. There is an apparent excess in the modeled [S II] emission, which may 
be remedied by reducing the assumed abundance of sulfur in the model (we used solar 
abundances throughout). Allowing for this additional free parameter is beyond the scope of 
this paper, however, as we are more interested in showing a consistency in the line profiles 
rather than their relative strengths. 

Figure 2 shows the radial dependence of the gas temperature and the ionization fraction 
of hydrogen in the best fit model. From this we see that most of the hydrogen is ionized 
near 1 pc from the center. So significant free-free absorption is expected in the region. 

Thermal equilibrium occurs when the heating and cooling rates balance; stable 
equilibria occur only where the slope of the cooling function is positive (e.g., Shore 1992). 
Figure 3 illustrates the heating and cooling functions for the gas at radii r = 10 17 cm (top) 
and r = 10 20 cm (bottom). The heating rate includes both gravitational dissipation and 
radiative effects, and the cooling rate must be determined including the dependency on 
the ionization state of the gas. Cloudy calculates both of these rates in a self-consistent 
manner. At r = 10 17 cm, it is clear that only one stable solution exists, T ~ 10 5 K (compare 
with Fig. 2), so there is no ambiguity in determining the equilibrium temperature. At large 
radii, however, it is apparent that two valid physical solutions exist. Because our picture 
has the gas being heated from a low temperature as it accretes, we assign the smaller value 
in all such cases of uncertainty. As shown in figure 2, the gas temperature jumps from 
~2x 10 3 K to ~ 10 4 K near 4 x 10 19 cm when the gas moves from the lower temperature 
stable branch to the higher temperature one. This is interesting because neutral atomic 
hydrogen is concentrated at that radius and this may therefore explain the detected neutral 
hydrogen radio line absorption. 

We have shown that our model can reproduce the observed broad-line spectrum; 
we now turn to the question concerning the observed absorption features. The free-free 
absorption optical depth at frequency v is given by (Walker et al. 2000): 



where all the quantities are given in cgs units. The line-center H I A21 cm optical depth, 
corrected for stimulated emission, is given by (Rohlfs, Wilson, & Huettemeister 2000): 




(9.8 x 10~ 3 ) n\ IT-^v- 2 



17.7 + In (T L5 v~ v ) 



(27) 




1.45 x 10~ 15 cm 2 Hz K 




(28) 
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where N H is the neutral hydrogen column depth and 




2kT 




(29) 



is the Doppler line width (including both turbulent and thermal broadening). Note that for 
r > r c , the turbulent velocity of the gas is always much greater than the thermal velocity. 
Figure 4 shows the radial dependence of both of these optical depths for our best-fit model, 
under the assumption that the line of sight is perpendicular to the disk. The free-free 
absorption is strong enough to account for the observed radio shadowing. However, the 
neutral hydrogen absorption seems deficient. This is not surprising given that we have 
neglected the density contrast of the supersonic disk in the calculation. Observations 
of other radio galaxies have revealed that the gas is clumpy at the same spatial scale 
(e.g., Harms et al. 1994; Haschick, Baan & Peng 1994; Kartje, Konigal & Elitzur 1999). 
Introducing density contrast will not only change the temperature profile of the gas, but also 
increase the neutral hydrogen column density. So we would expect significant enhancement 
of radio line absorption in a detailed modeling that incorporates these additional effects. 

Figure 5 illustrates how the ionizing flux modifies the line spectrum. With an increase 
in Li on , the heating and ionization parameters of the gas increase as well. This has the effect 
of moving the emission zone to larger radii, particularly for the low ionization [N II] and 
[S II] species. As a result, the line profiles become narrower for increasing L ion . A reduction 
of Lion from the best fit value does not have such a drastic impact, however. This is because 
the inner radius reaches an asymptotic limit in the low luminosity case, determined solely 
by balancing the gravitational heating and line cooling rates via Equation (15). In this 
limit, the additional heating due to the radiation field is negligible. It is interesting to note 
that our best fit lies on the transition to this limit. 

The dependence of the line shape on the stress parameter a is demonstrated in Fig. 6. 
A larger a increases the accretion rate (for a fixed number density at r a ), causing additional 
heating in the gas (via Eq. 13). This moves the emission zone to larger radii, resulting in 
narrower line profiles. Decreasing a has the opposite effect. 

We have used the observed radio free-free absorption to fix our value of n Q in the best 
fit model. Figure 7 illustrates the effect of varying n . The density of the gas is important 
in many aspects, including a determination of the ionization parameter, as well as the 
gravitational heating and line cooling rates. It also has a strong impact on the fluxes of 
the [SII] and [Nil] forbidden lines, which are suppressed relative to Ha for densities greater 
than the critical density, at which the radiative and collisional transition rates balance. 



Because the line emission broadening is due mainly to the Keplerian motion of the 
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gas, the line width is proportional to t> fc sini, where Vk is the Keplerian velocity in the disk. 
The projection effect, characterized by the inclination angle i, is demonstrated in Fig. 8. 
As expected, model spectra associated with small inclination angles have very narrow line 
profiles. Note, however, that even when % — the lines have non-zero width. This is the 
effect of broadening due to the turbulent velocity v t = \pf v k- 

Finally, the turbulent parameter / has a similar effect on the accretion rate as the 
stress parameter a. However, by increasing /, we not only shift the broad emission line 
region toward large radii, but also increase the total line flux, as the disk becomes thicker 
(thereby increasing the size of the emission region). Figure 9 shows how the line spectrum 
changes with /. 

5. Discussion 

The accretion flow within rt is not well constrained. The fact that the inferred accretion 
rate from this supersonic disk model gives a radiative efficiency of ~ 2 x 10~ 4 suggests that 
advection may dominate in the inner region (Narayan & Yi, 1994). However, if we associate 
the observed X-ray emission with such an advection dominated accretion flow (ADAF) 
at small radii, we find that the observed X-ray spectrum with a spectral index of ~ 0.4 
(Rinn & Sambruna 2001; Matsumoto et al 2001) is somewhat softer than that of the model 
prediction (Narayan et al. 1998). An ADAF with strong winds may be able to produce a 
softer X-ray spectrum to account for the high energy emission from the AGN in NGC 4261 
(Quataert & Narayan 1999). And this outflow may be associated with the jet formation 
process in the core. 

It is also interesting to note that NGC 6251 is quite similar to NGC 4261 in that 
both of them have a hundred-pc scale dusty disk, large scale symmetric radio jets and an 
inverted radio spectrum at the core (Melia, Liu & Fatuzzo 2002). Compared with NGC 
4261, NGC 6251 is about 3 times further away from us, its X-ray luminosity is 30 times 
stronger, and its broad Ha line flux is about 100 times bigger. HST observations suggest 
that there is a supermassive black hole of 6 x 10 8 M , which is very close to the black 
hole mass in NGC 4261. If a similar accretion flow exists in NGC 6251, Equation (21) 
suggests that the electron number density must increase by a factor of ~ 10 to produce the 
strong Ha emission. Combining this with the fact that the disk has an inclination angle 
of ~ 32° (Ferrarese & Ford 1999) and that the AGN is more active, we would expect that 
the observed broad emission line may also be fitted with our model. We also would expect 
strong radio line absorption in the direction of the counterjet with an H I column density of 
~ 10 23 cm" 2 . 
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In an earlier paper we introduced a hot expanding gas model for the mult i- wavelength 
spectrum of the AGN in NGC 6251 (Melia et al. 2002). According to that picture, a 
standard a-disk with an accretion rate of ~ 4 x 10 22 g s _1 exists near the black hole. How 
such a fossil disk can form after the supersonic disk terminates is still an open question. 
However, the implied strong evaporation of the accretion flow below r t suggests another 
origin for the hot unbound gas proposed in that model. Further investigation is clearly 
warranted, and an application of the supersonic disk model to other radio galaxies may also 
prove to be invaluable in developing a unified theory for supermassive black hole accretion. 
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Fig. 1. — Observed (points with error bars) and best-fit modeled (solid curve) broad optical 
line spectrum of NGC 4261. The best-fit free parameters are L ion = 4 x 10 40 erg/s and 
a = 0.4, resulting in xt — 11 iy — 250 degrees of freedom). The other parameters are fixed 
in the fit: n Q = 30000 cm" 3 , / = 0.045, and % = 64° (see text for details). 
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Fig. 2. — The electron temperature (left) and ionization fraction (right) profiles of the disk. 
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Fig. 3. — Cooling (solid curves) and heating (dashed curves) rates in the disk gas at 
r = 10 17 cm (top) and r = 10 20 cm (bottom). Note that the equilibrium temperature 
occurs where the heating and cooling rates balance. Stable equilibria occur only where the 
slope of the cooling function is positive. See text for details. 
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Fig. 4. — The free-free optical depth at 5 GHz (solid curve) and the line-center optical 
depth TQ lcm of H I A21 cm (dashed curve), calculated for the supersonic disk viewed face-on. 
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Fig. 5. — The dependence of the broad-line spectrum on the ionizing luminosity Li on - 
Shown in the figure are the modeled spectra for L ion = 4 x 10 39 erg/s (dot-dashed curve), 
Lion = 4 x 10 40 erg/s (solid curve), and L ion = 4 x 10 41 erg/s (dashed curve). The other 
parameters are fixed at: a = 0.4, n = 30000 cm~ 3 , / = 0.045, and i = 64°. 
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Fig. 6. — The dependence of the broad-line spectrum on the stress parameter a. Shown in 
the figure are the modeled spectra for a = 0.04 (dot-dashed curve), a = 0.4 (solid curve), 
and a — 1 (dashed curve). The other parameters are fixed at: L ion = 4 x 10 40 erg/s, 
n = 30000 cm" 3 , / = 0.045, and % = 64°. 
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Fig. 7. — The dependence of the broad-line spectrum on the number density n Q of the disk 
at a distance of 10 18 cm from the black hole. Shown in the figure are the modeled spectra for 
n = 3000 cm -3 (dot-dashed curve), n Q = 30, 000 cm -3 (solid curve), and n Q = 300, 000 cm -3 
(dashed curve). The other parameters are fixed at: L ion = 4 x 10 40 erg/s, a = 0.4, / = 0.045, 
and % = 64°. 
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Fig. 8. — The dependence of the broad-line spectrum on the inclination angle of the disk 
i. Shown in the figure are the modeled spectra for % = 0° (dotted curve), % = 30° (dashed 
curve), % = 64° (solid curve), and % = 75° (dot-dashed curve). The other parameters are 
fixed at: L ion = 4 x 10 40 erg/s, a = 0.4, n Q = 30000 cm" 3 , and / = 0.045. 
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Fig. 9. — The dependence of the broad-line spectrum on the turbulence parameter /. 
Shown in the figure are the modeled spectra for / = 0.011 (dot-dashed curve), / = 0.045 
(solid curve), and / = 0.180 (dashed curve). The other parameters are fixed at: L lon = 
4 x 10 40 erg/s, a = 0.4, n a = 30000 cm" 3 , and % = 64°. 



